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\0 ' Abstract 

o ■ 

OO ' The structure of the QCD gluonic cascade in configuration space is investigated. The 

exphcit form of the inclusive single particle density in configuration space transverse coor- 
dinates is derived in the double logarithmic approximation (DLA) of QCD. The possible 
simplification of the multiparton density matrix formalism for DLA approach is found and 



Q I discussed. 
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1 Introduction 

The aim of this talk is to present the recent results of the investigation |l|] of the 
structure of the QCD gluonic cascade in configuration space. The interest of the 
configuration space structure of a hadronic source has appeared primarily in intensity 
interferometry . The technique was developed originally to estimate the dimension 
of distant astronomical objects. Since that time it has seen widespread application in 
subatomic physics, in particular in analysis of elementary particle collisions [0. The 
standard HBT procedure involves introducing an ansatz describing the geometry 
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of the particle source, on the basis of a physical model. Then it investigates its 
multiparticle characteristcs in momentum space. Many models of hadron production 
based on some phenomenological or theoretical constraints have been considered; 
however, the question as to what is the configuration space structure of the QCD 
cascade when derived explicitly from the fundamental theory has not been addressed 
so far. 

Recently, several groups have analysed in great detail the multiparton distri- 
butions in the QCD gluonic cascades |Q. The results of their investigations show 
that perturbative QCD p provides a powerful framework not only for the descrip- 
tion of hard quark and gluon jets but also of much softer multiparticle phenomena. 
Although not understood theoretically, the hypothesis of parton-hadron duality 
provides an apparently successful hnk between theoretical parton distributions and 
observed particle spectra. This prescription was extensively tested in single particle 
spectra ( and total multiplicities) and found to be in a good agreement with the 
available data (see e. g. 0). Recently, there have appeared indications that it may 
also work for multiparticle correlations 0. These unquestionable successes invite 
one to study further consequences of the theory for processes of particle produc- 
tion. At this point we would like to notice that, if one wants to exploit fully the 
quantum-mechanical aspects of the QCD cascade, it is necessary to study at first 
the multiparton density matrix. The multiparticle distributions calculated so far 
give only diagonal terms of the density matrix and thus represent a rather restricted 
( although very important ) part of the information available from the theory. 

It is perhaps important to stress that in contrast to what is usually believed, 
the interest in studying the multiparticle density matrix is not purely academic. As 
suggested in 0, the density matrix allows one to obtain the multiparticle Wigner 
functions and consequently gives information about the space-time structure of the 
system. This allows one to make predictions about the shape and range of the HBT 
interference with clear experimental consequences. 

In my previous paper |T^ I investigated the multiparticle density matrix (DM) 
of the gluonic cascade produced in e~^e~ collision in the framework of double loga- 
rithmic (DLA) approximation [Q, |Tl|] of QCD. I proposed a generating functional to 
obtain integral equations for the multiparticle density matrix in the quasi- diagonal 
limit, i. e. if the energies and emission angles of particles are close to each other. 

Here I will be presenting the technique for extracting physical information from 
the density matrix approach, deriving the explicit form of inclusive single particle 
density in configuration space. For the sake of simplicity, I will be restricting myself 
to a discussion of its dependence on the transverse space coordinate xt- In solving 
the problem, first I will discuss the derivation of the inclusive single particle density 
matrix d^p{k',k). I will prove that the terms in d^p{k',k) which do not vanish for 
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k = k' are leading in the double logaritmic perturbative expansion when 05 ^ 0, 
momentum of quark ( antiquark ) P generating the gluonic cascade is large : 
P ^ 00 and the double logarithmic contributions of the form In P = const 
dominate. Taking into account these leading terms will allow me to obtain the 
explicit analytic form of the density matrix, and consequently the explicit form of 
the single particle density in configuration space. Finally, I will be analysing the 
physical properties of inclusive single particle density, and I hope to find them in 
agreement with intuitive physical expectations. 

2 Density matrix formalism 
2.1 Definition of the density matrix 

In this paper I will concentrate on the single particle density in configuration space. 
I will be discussing it in relation to the QCD-parton cascade within the framework 
of double logarithmic approximation (DLA), using the density matrix formalism 
presented in [|rU]. Hence in this section I will recall the definition of the density 
matrix for a particle production process, and show the relation between the density 
matrix and particle density in configuration space. Our assumption is that all the 
produced particles are real, i.e. they are on a mass shell. If the production of 
m particles can be realized in different ways represented by a sample of Feynman 
diagrams, then the exclusive m-particle density matrix equals the product of total 
production amplitude S{ki, . ., km) and its complex conjugate *S'*(/b^, . ., /c^) as '. 

d (k^, . ., k^', ki, . ., km) = S {ki, . ., km)S{ki, . ., km)] (1) 

where the total amplitude S{ki, . ., km) is the sum of all contributions from 
Feynman diagrams (D) multiplied by phase space factors (2co'fc)"^^^ : 

m 

Sik,, . ., U = E S(D)ih, ..,km) l[{2uJkr'/^. (2) 

D i=l 

For inclusive analysis one constructs the m-particle density matrix as a series of 
integrated n-particle exclusive densities in the form : 

r/^fk' k' ■ ki k — 
°° 1 /• 

E 7 ~ [d^m+i.^n c^^^lk'i, . ., k^, km+i, . ., k„; ki, . ., k„); (3) 

n=m TTl). J 

where [dk]i,,,j = d^ki. .d^kj. This construction scheme implies, that the diagonal 
elements of the density matrix are equal to particle densities in momentum space. 
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There is also an obvious relation between the density matrix and particle density 
in real space. Remembering, that the space-time multiparticle amplitude is the 
on-mass-shell Fourier transform of the momentum amplitude, i.e. 

where denotes the energy of the ith produced particle, for the multiparticle 
density p^^/*"(xi, . we get the relation : 

p^^/^"(a:i,..,a;J = ^^^ / [rfA;]i...^[rffc']i.,.^ . ., ^ fci, . ., U (5) 

where the factor ^^^-^^3^ was introduced to get the proper normahzation : 

2.2 Parton cascade in the double logarithmic approximation 

The double logarithmic approach in momentum space gives a good quali- 

tative description of the structure of the gluonic cascade. It accounts only for the 
leading DL contributions to the multiparticle cross section. Although emitted soft 
gluons violate the energy and momentum conservation rules, however, at high ener- 
gies the approximation reproduces quite well the selfsimilar structure of the gluon 
radiation. Let us consider in the framework of DLA the gluonic cascade generated in 
e+e~ collision. Multiparticle exclusive amplitude S'e^. .e„(/ci, • •, /cm), describing the 
production of m gluons : 

5e,..em(fci,..,U = (-ire-^nMp,(fc,)e-^; (7) 

i=l 

is a product of the m emission factors : 

MpXK) = gs^^^^QpAki)Gp- (8) 

where: 

gs = A/47raI, 

n is the number of gluons emitted of quark (antiquark), 
ki = (tUj, kj) denotes the 4-momentum of the ith soft gluon. 
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d = Cj , j = 0, . . . , 3 describes its polarization, 

Pi is the 4-momentum of the parent of the ith gluon, 

Gp. denotes the color factor for the given vertex of the tree diagram D, 



and can be conveniently represented as a tree diagram (see Fig. 1). Gluon emissions 
are not independent. Radiated particles have the memory of their parton parent, 
and of the previous parton splitting off its parent line. This dependence restricts the 
phase space of the produced gluon, and it is included in the form of a generalized 
step function 0p.(/ci) : 

Opik) := {k^ = co< P°, ^^kP < 0, u9kp > Qo}; (9) 

where P is the momentum of the parent of a given parton k, 9 denotes the emission 
angle of the previous parton splitting on the P— hne, and Qq is a cut-off parameter. 

■w(P) 

Virtual corrections appear as a radiation Sudakov factor e 2 ^ where : 

w{P) = J d'k{AUk)Ap{k))^,) (10) 

denotes the total probability of emission of a gluon from the parent P, averaged over 
physical transverse polarizations e^,e^. Ap{k) is given by : 

, , , Mp(k) , , 

It should also be emphasized that produced gluons are real (on-mass-shell) particles, 
so the energy a;^ of a gluon of momentum k can be approximated as : 

cufc =1 k I . (12) 

Summing of the color factors G over the color indices gives the result : 

^pAp. - j Cv dla Pi ^ P; ^^^^ 

where P denotes the 4-momentum of the quark (antiquark) which initializes the 
gluonic cascade. 

In DLA different tree diagrams come from different non-overlapping kinematic re- 
gions, and do not interfere. Therefore, to calculate exclusive and inclusive multigluon 
densities it is enough to sum up incoherently the squares of amplitudes (0). Hence 
one obtains the exclusive density pp (ki, . .,km) in the form : 

m 

p-(ki,..,kj =e-(^) ^n(^K(^^)^P.(^^))(e.)e-"'^''^ (14) 

D i=l 



5 



parametrized by the momentum P of the quark (antiquark) which initiahzes the 
cascade. Multigluon inclusive density pp (ki, . .,km) follows from as : 

°° 1 r 

Ppiki, . ., k^) = 1—^ — \\ / Pp (ki, • k„). (15) 

n=m, ^J- 

Introducing the method of the generating functional (GF) (see and references 
therein) allows us to perform the summation over diagrams in ([T^ and (^) in 
a very convenient way. While constructing the generating functional Zp[u\, one 
applies explicitly the self similarity property of the gluonic cascade. As a final result 
one obtains the recursive master equation in the form : 

Zp[u] = e-"'^^) exp{ f d^k {A*p{k)Ap{k))^e) u{k)Zk[u]). (16) 



It can be proved that equation (|T^) reproduces contributions of all tree diagrams 
D, and allows us to express multigluon densities pp(ki, . .,km) and pp (ki, . .,km) 
as : 

p-(ki, . ., k^) = ^Zp ||,=o}, (17) 

OUi . . . OUm 

p-(ki, . ., k^) = ^Zp (18) 

OUi . . . OUm 

where Ui denotes the probing function u{ki) and the functional derivative J- denotes 



respectively. 



5u{ki) 

We would like to emphasize the simplicity of the GF approach. Especially for 
inclusive densities the method allows us to skip the complicated summation proce- 
dure, and to express the required distribution in a simple, compact form (for details 
see 1112 



2.3 Density matrix in the DLA formalism 

The DLA formalism in momentum space gives a good description of the structure of 
the gluonic cascade. The GF scheme suggests clearly how to construct multigluon 
densities, and allows us to investigate their properties in a simple way. Since the 
method works so well for multiparticle distributions, one can expect to apply it 
successfully for other multiparticle observables. Below we briefly discuss the calcu- 
lation of the multiparticle density matrix in the framework of double logarithmic 



approximation, recalling our main results from Ref. [jTO . 

First we derive the general expression for the exclusive and inclusive density 
matrices d^p(k[, . ., k^; ki, . ., k^) and rfp (k'^, . ., k^; ki, . ., k^)- The task looks quite 
complicated because in this case the interference between different diagrams in (Q), 
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generally does not vanish. Let us define two functionals, Zp[u] and Zp,[w], which 
generate the sum of all tree amplitudes and the sum of their complex conjugates 
respectively : 

Zp[u] = e-'^^^^/^expiJ d^kAp{k)u{k)Zk[u]), 

Z*p,[s] = e-""^^'^'^ exp{j d^kA*p,{k)s{k)Zl[s]). (19) 

The multigluon density matrices can be then expressed as (see Appendix A) : 
dp • - I ^^m) ki, . ., km) = 



5Sii . . . 6Sm' Sui ■ ■ ■ 



dp 0^1, ■ •) k^; ki, . ., km) 



-Zp[u]Zp,[s] |{„=5=o},p=P' 5 (20) 



= 1 X 1 X Zp[u]Z*p,[s] I r A 01 p^p, . (21) 

6Si' . . . 5Sm' 6Ui. . . 6Um ^ L J 4n g^,s u>,i- 

Eq. ([21| ) is, in fact, a complicated integral equation. The difficulty of the diagram 



summation does not disappear there. It is only hidden in the compact form of (pij). 
However, we remember that in the above formula we have taken into account inter- 
ference between all different graphs D and D'. And detailed analysis gives the result 
that in DLA not all the diagrams mix up: one can distinguish some interference 
classes. However, at the general level we did not succeed in formulation of such a 
GF which would take this fact into account. 

Nevertheless, we do not need the most general form of the density matrix. We 
are interested in its behaviour when the differences of momenta 
I ki — k[ I , • • • , I f^m — k'^ I are small, since we expect that large momentum 
differences will not contribute to Fourier transforms [|l^]. It can be shown that 
within this limit interferences between different diagrams vanish, and one sums up 
only "squared" contributions from identical graphs. In fact, from the analysis of 
the diagrams contributing to single particle density matrix (ip(fc^,fci), it has been 
possible to prove (see Appendix B) that interference between different diagrams 
appears only if either energies or emission angles of produced particles are strongly 
ordered : ui ^ LUy {lji <^ uji>) or 9ip ^ Oyp {Oip ^ Oyp). This statement can 
be generalized for any m-particle density matrix. If we have m particles, and (m-1) 
ones among them are "close" to each other, e.g. ki = k[, . . . , k^-i = ^m-i' then 
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the interference of the different diagrams will take place only if either the energies 
or the angles of and k'^ are strongly ordered. 

Hence in our approximation we exclude mixing up different diagrams, and sum 
up only the squared contributions from identical ones. The exclusive and inclusive 
density matrices then take the simpler form : 

'^p^(ki, • k^; ki, . ., km) = 

m 

= E n(4^fc^^^J-'/' {Se...e^{k[, . ., 0^e...e„(fcl, . ., U)(e...e„), (22) 
D i=l 

dp (ki, • k^; ki, . ., km) = 

oo -1 „ m n 

= E 7T— TV E/M^]™+i...nn(4^M^^^)"'^' n (4^..^..)-^/^ 

n=m \ I'D i=l j=m+l 

X {Sei. .e„{ki, . . . , k^, kfn+1, ■ ■ ■ , kn) Se^. .e„{kl, ■ ., ^n) ) (ei. .e„) • (23) 

The summation in (|23|), (|23D over D can be easily performed using the generating 
functional which reproduces only contributions of identical tree diagrams. In ||lO| I 
proposed a master equation for such a generating functional (GF) Zp'p[u{k', k)]. It 
accounts only for the contributions of "diagonal" diagrams (for proof see Appendix 
C), and generates the exclusive and inclusive density matrix as : 

(ip (k^, . ., k^; ki, . ., km) = Zp/p \ fu=o},p=p', (24) 

dUm',m- -OUif^i 

dp (k^, . ., k„; ki, . ., km) = t j ^p'p \{u=5'i{i'-i)},p=p' ■ (25) 

The explicit form of its master equation, and the other formulae needed for further 
calculations are presented below. The GF reads : 

oo 1 

Zp,p[n] = e-^(^''^)^- / [dk']^,„n[dk],.„^uik[,k^)..uik',,kn) 

x{A*p,{k'Mpiki))ie,)--{A*p,ik'JApik^))^,^) 

X^fcJfciM- ■Zk'„kju] Pl',..,n';l,..,n; (26) 

where the function Pi',. .,n';i,. .,n provides the requested parallel angular ordering ( see 
Fig. 5) for n particles in the form : 
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Pl',..,n';l,..,n = ©(^'fe' p' > . . > 6k' P') Q{6ki^p > . . > Ok,„p)] (27) 

(il,. .,in)G-Pe'"'Ti(l,- -j™) 

the product of the single particle amplitudes ( A*p,{k')Ap{k))(^e) averaged over gluon 
polarizations (see Appendix D) reads : 

{A*p,{k')Ap{k))(^e) = Ap,p{k',k) = 

= ^ G^Gp ^^^_0p,(A:')ep(A:); (28) 
ZTT VAuj'-^uj'^ ypk"P'k' 

and the radiation factor is given by : 

M,'(P',P, = ^m±^. (29) 

The other notations are the same as in section 2. 2. 

Generating functional ( p6D describes the sequence of "emissions" : 
(Ap,(A;')y4p(A;))(e) of two particles fc, k' from two parents P, P' . If the profile function 
u{l',l) is equal to and 6^{l' — I) respectively, then the functional Zp'p takes the 
form : 

Zp,p[n = 0]= e-^(^''^), 
Zpp[u = 6%l' - 1)] = 1. (30) 

For the profile function of the form u{l', I) = v{l) 6^{l' — I) and P = P' our Zpip[u\ 
reduces to functional (0), as expected : 

Zpp[u{l\ I) = v{l) S\l' - I)] = Zp[v{l)]. (31) 
2.4 Single particle density matrix. Leading terms 



From relations (pSf), ( [26| ) one obtains a simple integral equation for the inclusive 
single particle density matrix : 

d'^ik'-k) = j (fsApp{s,s)d'^\k';k) + 

+ fp{k',k)App{k\k)Zy,k[u = 5\l' -I)]- (32) 

where fp{k',k) equals 

fp{k',k) = e-^sC.\ir^''^-'''''-^l (33) 
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Introducing the notation : 

gp{k', k) = fp{k', k) App{k', k) Zk',k[a = - /)]; (34) 
we may write the symbohcal solution d^p{k' , k) of (p^ in the form : 



oo „ 

(f^{k',k) = J2 d^si. .d^SnApp{si,Si). .A,^_^s„_^{sn,Sn) gsAk',k). (35) 

ri=0 



We present Eq. ( |35D as a final result obtained in Ref. [|TOl- The exact solution 



of d^p{k',k) and its detailed properties will be discussed below. We may replace 
App{s, s) in ( p5D with its explicit form taken from Eq. (p8|): 

d^p-{k',k) = £(2^/— //§^ep(.o../— 

X9sAk',k); (36) 

where / dQsi denotes the angular integration over the direction of vector Si, Bjj = 
Qsi,sj and b = g'sCp- Factor gs{k' , k) defined by expression (P^), reads : 

X e "30 Qo ' e Qo Qo ' 

°° , f 1 111 

X YS^b'r J [dk],...r. i^^^e^^- ■ 2^e^j^,^ ^- 

X 6(^1 > . . > efcn)0(^fc'l > • • > Ok'n) 

X 6^,(1). .efc,(n)0fe(l)..0fc(n); (37) 

where V = ggCy- Rewriting its phase space restrictions in terms of Bessel functions 
(see Appendix E) one arrives at : 

X e ' «o ' e «o + <5o ^ 
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" / roo rmin(k,k') rjh. 



n=0 j=l 



X d9kiJo{xi9ki) dOk'iM^i^k'i)) ■ (38) 



m,aa;(3 ) k 

"k.i-iik'i-i) 



The term which dominates in ( p7D in the double logarithmic perturbative ex- 

Qo 



pansion : b ^ 0, P ^ oo, 61n^^ = const, has the explicit form (for details see 



Appendix F ) : 



9^s\k', k) = ^ ^s{k). (39) 



One may check (see Appendix F ) that this term iterated in (|36D produces the DLA 
leading contribution to the density matrix d^p{k',k). Namely, if we introduce the 
notation : 



n=0 J Si J ZTTU^p^ J Sn J iTlU^^i^^ 

X9i'Jik',k); (40) 
then there is a relation : 

(r^{k',k) < d!p'^^\k',ky, (41) 

and for k = k' in quasi- diagonal limit : 

(fp{k',k) ^ df^^\k\k). (42) 

Both ( ^iD and ( ^21) hold in the perturbative limit, as well. Hence we may write 
finally that : 

°° - dsi f dQsi ^ f dSn f dQ^n ^ / ^ 

c„ T [Sr. 



^ J Si J J Sn J 



2 ^Sn-iK-^nJ 



PI J Sn J ^nOn-i 

xgW{k\k). (43) 

The above result simplifies significantly the single particle density matrix approach 
p!0| . We apply it in section 4. 2. so as to improve calculation of the single particle 
density in mixed coordinates. 
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3 DLA in configuration space 



3.1 Fourier transform in transverse coordinates 

The DLA soft parton cascade starts from the initial parton of momentum P. Mo- 
menta of all particles produced from parton P refer to the P direction: they depend 
on the transverse and longitudinal momenta kx and kL (see e. g. 0, [T^) taken 
with respect to the P axis. 

This dependence influences the structure of the cascade in configuration space. 
For the single particle amplitude Eq. (^) takes the form : 

Sp{x) = J d^ke^'^'-'^^'Spik) = I dkLd^kre'^^'-'^^'Spiky, (44) 

where 5'p(k) {Sp{x) ) denotes the amplitude to produce a single particle with the 
momentum k ( at the space-time coordinate x ) from the initial particle of momen- 
tum P, and the 3-dimensional product kx reads : 

kx = fc^Xi + kx xt; (45) 

where indices L and T denote longitudinal and transverse components of momenta 
and coordinates with respect to the P axis. 

Since DLA amplitude Sp = Sp{uj,9pk) depends explicitly on u; =| k |= 
(p!^) and 9pk, one should express it as a function of ki and ky. To proceed, let us 
recall that the simple form of Sp = Sp{k, 9pk) has been obtained from the original 
QCD expression using the approximation of small emission angles <^ 1 [Q. In 
this approximation one retains only the leading term in 9pk. Hence, restricting 
ourselves to the leading terms in 9pk-, we interprete product k9pk as the transverse 
momentum fc-p and k as the longitudinal momentum k^. The correctness of the 
substitution k ^ k^ for e*'^* in (^) is not so evident, since it should be done in the 
argument of Fourier transform. However, if we restrict ourselves to finite, small time 
intervals, i. e. t -C {P9'^)~^, it also works with a high level of accuracy. To see this, 
let us analyse the dependence of /c on /c^ and kx- Following (0), it can be expressed 
as : 

{uj=)k = ^Jkl + k^. (46) 

For intrajet cascades there is a relation kj- <C ki. Restricting ourselves to the 
leading term k ^ k^, we ignore the contributions of the order of -r^. In DLA the 

2 

ratio TT has a limiting maximum value 9"^ denoting the emission angle of the previous 
parton splitting off the P-line (compare @)- Maximal energy carried by the gluon 
equals P. Since that the approximation k ^ ki implies, that one neglects terms of 
the order of P9'^. 
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Using the leading term approximation, Eq. (|4^) can be rewritten as the product 
of 2 separate Fourier transforms: 1-dimensional FT of {ki, t — xi,), and 2-dimensional 
FT of (kx, — Xt) in the form : 

Sp{xL,^T,t) = J dkL(^^^^*-''^^ J d^kre-'^^'^^Spik ^ fci.kx). (47) 

Eq. (^Vp gives us a good tool to investigate the space-time structure of DLA. We 
may now concentrate on particle distribution in transverse coordinates, which is 
physically the most interesting case. For the sake of simplicity, let us consider 
for if: = the single particle amplitude S'p(/c,xt) in mixed transverse space and 
longitudinal momentum coordinates [jl4| defined as : 

Sp{kL ^ k, XT, 0) = y d^kr e-'^^^'^SpikL ^ k, kx). (48) 

For the single particle exclusive distribution pp(fc, xt,0), defined as the amplitude 
square : 

P^^(A:,xt,0) =|5p(fc,XT,0) |2; (49) 
one obtains simple expression in the form : 

p^^'(fc,XT,0) = j^ld'kTd'k'Te~'^'^^-''''r^^^Spik,kT)SUk,K) 

= — ^ /d'A;TC^'A;^e-'('^^-''T)-:rci^-(fc,k'T;A:,kT). (50) 

Relation (l50|) holds also for inclusive single particle density (for proof see Appendix 
G): 

P*p"(A;,xt,0) = fd'kTd'k'^e-''^''^~''''r^''^d'P{k,k'T,-k,kT:y, (51) 

and can be easily generalized for any multiparticle distributions. In this paper we 
study the inclusive single particle density in mixed coordinates pp (fc, xt, 0) for the 
QCD gluonic cascade in DLA approximation. 



4 Inclusive single particle density in configuration 
space 

4.1 Physics 

Momentum and configuration space description of a particle source are related to 
each other by the Fourier transform. Therefore one expects to extract from the 
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observables in momentum space some qualitative information about their behaviour 
in configuration space. This observation apphes, of course, to single particle density. 
In the DLA approach the inclusive single particle density Pp (k) ||12| reads : 

- Im'^ ep(.); (52) 

where h = gsCp. For constant k the Pp (k) is concentrated around P direction (see 
Fig. 6)( skipping for now the 0p(/c) limitations ). We expect, that the single particle 
distribution pp(t, x) in configuration space : 

f/Pik, XT, 0) = J (fkT(fk'T e-^('^T-k^)''^ 4"(A;, k'^; k, kx); (53) 

will be concentrated around the P direction, as well. However, it has to obey the 
uncertainty principle. Since there are cut-offs of form (^ in density matrix, its 
Fourier transform will contain some oscillations, due to the restricted integration 
region. Moreover, relation (^ implies, that both densities Pp(k) and pp(/c, xx,0) 
have the same normalization, if integrated over (Pk and dk (Pxt respectively. 

Remembering the above remarks, we propose now a technique for deriving the 
explicit form of inclusive single particle density in configuration space. 



4.2 Single particle density 

We may calculate the single particle density pp (A;, xt,0) as a Fourier transform of 
DLA density matrix (|43|). Let us make the transverse Fourier transform of both 
sides of Eq. (H). We obtain : 



p^"(fc,XT(P),0) = I rf2A;r(P)t^'A;^(p)e-^('^^(^)-'^np))''^(^)4"(A;,k'T;A;,kT) 



DLA 

- (2vr) 



X 



^/ou\n fdsi f dn,-^ r^ / \ f dSn f rffi,„ 

X E(2&) / — / T^QpiSl). . — 7-^2 0..-i(^n) 



?i=0 



Si J 27i9j,-^ J Sn J 27rei_i^^ 



X 



^?«(/.,k'T(s);fc,kT(s)); (54) 
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where indices T(P),T(S) correspond to the reference frame with the axis placed 
along the P [S) direction. The term of the lowest order [n = 0) takes then the 
form : 

Pp ^^^( A;, xt(p),0) = 

= ^^©(^ <k<P){ L ddpMxTkdpk)f. (55) 

The numerical plot of ( p5D is presented in Fig. 7. There is a limiting maximum value 
for xt = 0, as expected, and for large xt the function has a power decrease with 
the best fit exponent x^^'^'^ . The result confirms our intuitive analysis from section 
4. 1. The term pp '^^■'(/c, xt(p), 0) is concentrated around the P direction, however 
stronger than p^p^^\k, kx) from Eq. (|5^). 

Derivation of terms of an arbitrary order in expansion (^) is more complicated. 
To proceed, let us analyse the last part of (^). It reads : 

/ / 9^1?" Qn-l (gn) g^n^ (k, kfr(s) ; k, kT(S) ) = 



= VirJ—l erUsn)jr^enik')eM; m 

where the index (i) refers to the vector Sj. Taking into account phase space restric- 
tion and dominating contribution of g— — (see Appendix E) expression (0) can be 
approximated as : 

^ / — / l^^^Qn-ASn)j^Qn{k')Qn{k) = 

= -n e„_i(fc')0n-l(A:) X 

2lTk'^ k On-l,k'(^n-l,k 



X I dXnXnMXnOkk') jo^ dOnkJoiXnOnk) ' dOnk' MXnOnk') ■ (57) 



'n-l,fc' 



k 



Scheme (0) can be repeated iteratively in (^3]). After some transformations one 
obtains finally (for details see Appendix E, H ) : 

p^p"(fc,XT(P),0)= (58) 
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— k 0(-^ < k < P) dOpk dOpk' — dipk'k ■hi.kOk'kXT) 



X/, 



86 ln(— ) <j / dxxJo{x9kk') jg^ daJo^xa) da'Jo{xa') \) 



\ 'k'XJo 



^kQ{^<k<P) ll dOpk ll d9pk' ^ r ^J^'f''^ M kOk'uXT) 
2-K 9 J9^ J9^ 2-K Jo A{epk,6pk',9k'k) 



x/o(. 86 ln(— ) { / da I da' — — — 

"-k'^J^ J9a 2nA{a,a',9, 




where 6k'k denotes the angle between momenta k', k, which in the reference frame 
with the z-axis placed along P direction takes the simple form : 



Ok'k = V oik' + Opk-'^ dpk'dpk cos{ipk'k)] (59) 

(pk'k denotes the azimuthal angle between kx and k'rp, and A(a, b, c) denotes the area 
of the triangle with sides a, b, c. The numerical plot of (|58|) is presented in Fig. 8 . 
There is a limiting maximum value for = 0, as expected, and for large xt the 
function has a power decrease with the best fit exponent x^'^'^^. The result confirms 
our intuitive analysis from section 4. 1. The term pp^^"* (A;, xt(p), 0) is concentrated 
around the P direction, stronger than pp^^''(fc, kx) from Eq. (|5^) . The function 
oscillates around its power-law profile. As already mentioned in section 4. 1. , this 
effect is due to the sharp cut-offs of form (P) which restrict the integration space. 



4.3 Properties 

Let us analyse the properties of ( pHD in detail. Comparing ( |58D with the single 
particle density in momentum space (|5^) , we get the same dependence on the energy 
k, as expected. The complicated emission term (see Fig. 9) : 

26 In 1^— j <yj dxxJo{x9kk') J^^^ daJo{xa) daJo{xa) 

= I d's^As{k,k'T;k,kT) (60) 

corresponds to the logarithm In ^ of (|5^). In fact, after integration of p^p{k, xt(p), 0) 
over d'^XT one obtains the same result as for (^) integrated over d'^kT (see Appendix 
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I). The density p'^p{k, xt(p), 0) is positively defined, as well. Using the Bessel function 
identities from Appendix H, Eq. (^) can be rewritten in the form : 

p'p"(fc, XT(P),0)= (61) 
^0 N)' k Jo Jo mi..^=-oo 

X < d9pk Jrm{Xi9pk). -JmA^nOpk) Jmi + ...+mAxTk9pk) [[ J^^ dttiJo^Xiai) > ; 



which implies positive definitness. For xt = pp (fc, xt(p), 0) reaches its maximal 
value, as expected. Diagramatically, formula (^) represents the chain of indepen- 
dent emission factors ( |BI1| ) transformed onto the transverse xt plane by the Bessel 
factor of primary emission from the parent P, namely Jo(| — ^'t \ ^t) (see Fig. 
10). For k = P expression (^) reduces to (^) with the dominance of the primary 
emission, as expected. 



5 Summary 

We considered the QCD parton cascade created in e~^e~ collision in the double 
logarithmic approximation. Using the density matrix (DM) formalism |jTO[, and 



restricting ourselves to the terms leading in the double logarithmic (DL) pertur- 
bative expansion for the quasi-diagonal limit k' = k, we derived the explicit form 
of the inclusive single particle density matrix d'^p{k',k) and single particle density 
Pp (fc, xt(p), 0) (see Fig. 8) in mixed coordinates {ki ^ k^^x)- The gluon density 
Pp( A;, xt(p), 0) fulfills important physical requirements such as positive definitness 
and proper normalization. It is concentrated around the P direction, and shows the 
power law profile for large xj-. Due to the cut-offs of the type (^ which restrict 
kinematic regions, the density oscillates around its power law profile. 

The above results give a positive outlook for the future. The simplified technique 
for calculating multiparticle observables in DLA for the constant as is ready. It may 
allow one to investigate the structure of the QCD cascade in a very comprehensive 
way: the exact form of a Wigner function already would give us clear experimental 
predictions. 

In order to improve the approach one needs to include the momentum dependence 
of as- The exact analysis of applicability of the quasi-diagonal approximation 
k' = k would be required, as well. There is also a problem as to what extent the 
kinematical constraints characteristic for a given QCD approach ( in DLA they are 
of the form (^ ) infiuence the results obtained in the approach. In other words. 
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there is a problem how to separate the effects of the kinematical restriction from 
the dynamical content of QCD space. The explanation would require performing an 
analysis similar to the presented above for LLA and MLLA schemes, as well. The 
MLLA approximation would be of special interest since it would incorporate the 
energy conservation into the cascade. 

6 Appendix A 

The functionals Zp[u] and Zp'[s] defined in (|19]) generate the sum of amplitudes 
and their complex conjugates over all possible tree diagrams respectively. It can be 
seen when one rewrites the master equation for Zp in the form of the diagram series 
(see Fig. 2). From that construction follows exclusive density matrix (|20D : 

dp{k[, . ., k'^; ki, . ., k^) = Z*p[s] \{s=o} t j — Zp[u] \{u=o} ■ (62) 

dSi'...dSm' dUi...dUfn 

Then the inclusive density matrix calculated from (|^) looks like this : 




and from the identity for the product of any two functionals F, F' : 

F[u]F'[w] ||„.^,..o>= d'k--j F[u] |{„.o} F'H |{..o}; (64) 

follows (grp. 

7 Appendix B 

In DLA there are four different tree graphs Ma, M^, Mc, M^, describing the produc- 
tion of two gluons. They are defined on non-overlapping kinematic regions 
(see Fig. 3 ). Emitted gluons are either angular (AO) or energy ordered (EO). 

Let us consider all the diagrams contributing to the single particle density matrix 
d^p{k'i,ki). From the (AO) and (EO) it follows that the interference between any 
two different graphs will appear only if either energies uj\ , u'l or emission angles 
9ip, Oyp of produced gluons are strongly ordered (Fig. 4). 
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This statement can be generalized for any m-particle density matrix by induction. 
If we have m particles, and (m-1) ones among them are "close" to each other, i. e. 
ki = k[, . . . , km-i — k'^-i, then the interference of the different diagrams will take 
place only if either energies or angles of km and k'^ are strongly ordered. 



8 Appendix C 

Z p' p[u{k' , k)] defined in (^) generates correct exclusive density matrices (0). The 
proof of this statement follows from the representation of the master equation for 
Zp'p as a diagram series (see Fig. 5). The series reproduces all squared contributions 
of identical tree diagrams, and excludes interference of the different ones (function 

Pl',...,n';l,...,n)- 

The explicit form of inclusive density (^) follows from formula (0). Substituting 
(P^) into (H) one obtains : 

= J2-( <fk<fk'5\k' -k)-^] ^- Zp.p[u] |{.=o}; 



(65) 



which represents the functional Zpip[u] expanded around "null" {u = 0} in the 
"point" u = 6'^{l' -I). 



9 Appendix D 

We want to average relation (pS]) over physical polarizations. The exact expression 
to be summed over polarizations looks like this : 



{k ■ P){k' ■ P')' 

where e, e' are polarizations of the same gluon identical in the limit k' = k. The 
gauge fixing we apply in the approach allows us to neglect contributions to ( pBD 
coming from the " nonphysical" polarizations e° and e^. Furthermore, it requests 
the time components of the physical polarizations e^,e^ to be equal to 0. Space 
components of e^, can be then constructed in the form : 

2 _ P X k 
^ ~ iPxkl' 
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and for e' respectively : 



e'2 



P' X k' 



P' X k 



e'l = , (68) 
e'2 X k' ^ ' 



Summing over these 2 polarizations and expanding scalar products in (|66D, one 
obtains finally : 

,i(,2 (k-P)(k'-P') uiuj' OptOpv' ^ ' 

The result can be easily confirmed for any two physical polarizations e^, placed 
on the plane e^e^ : 

= cosLf + sinip , 

= —simpe^ + cosipe^] (70) 
and for e'^, e'^ placed on the plane e'^e'^ with the phases ip' = Lp respectively. 



10 Appendix E 



To transform (|56|) into (|5^) we will apply the pole approximation method ([0 and 
references therein). The phase space restrictions 9„_i(n)0„(fc)0n(^') in (^), in 
particular the angular ordering (AO) 6n-i-n ^ On,k, Gn-i-n ^ Gn,k' make the the 
term ^ — ^ — be dominating in (|56D . Applying the Bessel function identity fl^ : 



dfln 1 f dOnk dOnk' ^nk ^nk' 



2n 2n J A{6nk,dnk',0kk') 

dxxJo{x6k'k) / dOnkdnk JoixOnk) I dOnk'^nk' Joi^dnk')', (71) 



where 9k'k denotes the relative angle between vectors k, k', and A{9nk,0nk',0kk' 
equals the area of triangle with sides 9k'k, (^k'n, ^kn, we may rewrite expression (^ 
in the form : 



2b 
27rfc3 



/ — / ^#^0n-lW^0n(fcO0n(fc) 
J Sn J ^7ry„_i „ t/nk'fnk 
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2b f dSn. 



2nk^ J Sn Jo 

1 



- n2 Qn-i{n)Qn{k')Qn{k). (72) 

"n-\,n 

Because of angular ordering the angle 9n-i^n practically does not change while in- 
tegrating over angles Onk,Onk', and can be successfully approximated by the angle 
6n-i,k {0n-i,k')- Applying all the integration restrictions explicitly, we arrive finally 
at the expression : 

2S^/!f/2^^-^-)^^^(^')«"(^) - 



^ -e„_i(A:')e„_i(fc) T'^^x (73) 



27rk^ 9n-l,k'0n-l,k Jk s 

X / dXnXnJo{XnOkk') ddnkJo\XnOnk) L„ dOnk' Jo\XnOnk')\ 







k k 



which after trivial integration over s„ gives identity 



11 Appendix F 

We shall find such terms in expansion (|36| ) which dominate in the double logarithmic 
perturbative limit of d^p{k',k), i. e. when b ^ (6 oc as) and P oo, so that 
61n^ ^ remains constant, and generate double logarithmic corrections to the cross 
section. First let us introduce the following notation : 



oo 



n=0 



P''" / WJ 2^?tQ-(^''--/ if y ^e,._.(a„)9.„(.',.). (75) 



where coefficients a„ describe the nth order iteration of gs{k', k) (|36D in the form : 

dsi f dQg-^^ ^ / \ f dSn f dQ 
W 

The term oq of the series ( |75D is equal to gp{k', k). Expanding it in the powers of 
coupling constant 6, one obtains the term of the first (lowest) order of b in the form : 

° ^ ' ^ 27r v^Fps QpyQpk ^ ' ^ ^ 
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Denoting the other terms by the symbol a^^^^ = g^p^\k\ k), the term from the 



series (^5]) can be rewritten as : 

ao = 4'^ + (77) 

We have checked that a^^^ dominates in (^). To see this, one should rewrite ao as : 

ao = gp{k',k) = ai'^ + aS,>'^ (78) 
= a^^^ f{k\k-P)- (79) 

where the exact form of /(A;', k; P) follows from (0). For 6 > and P finite the 
f{k',k;P) vs I k — k' I is a gauss-like function with the maximum equal to 1 for 

PT 

k = k' and a non-zero minimal value. In the DL perturbative limit : f{k', k; P) 1. 
Consequently, for any k, k', P : 

ao < a^o\ (80) 
ao PT 



a. 



(1) 



1. (81) 







Since ao^"* dominates in ao, we expect that the iterations of ao^"* will generate the 
leading contributions to the dp (fc', k). Let us introduce the notation : 

oo 

d'^'^'\k\k) = ;^aW, (82) 

n=0 

oo 

<f-^^>'\k',k) = 5:a(>i); (83) 

n=0 

where a!^^^ and a!^^ represent the result of nth iteration (^) of ao"^^"* and ao^'' re- 
spectively. From ( pOD immediately follows the relation : 

an = a« + a(>i) < a^. (84) 

Since a„ > 0, for the density matrix one obtains : 

< (f^{k',k) = df^^\k\k) + df^''^\k\k) < d'^'^^\k',k). (85) 

In the limit | k — k' | ^ the density matrix d^p{k', k) d^p'^^\k', k), and there- 
fore : 

"1. (86) 
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For other ( larger ) values of | k — k' | the contribution of dp {k',k) generally 
need not to dominate. However, taking into account the fact that the density ma- 
trix approach works only for the quasi- diagonal limit, we may identify the quasi- 
diagonal region with the region of the dominance of (fp'^^\k', k). Furthermore, since 
(fp'^''^\k,k) = 0, the (fp'^^\k',k) generates all DL corrections to the cross section. 
Hence we finally arrive at the relation (^). 
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The particle density pp(fc,xr,0) can be represented as : 
°° 1 r 

Pp(fc,Xr,0) = ^ — — -— / dk2d'^X2T ■ ■ ■dknd^XnTPp{k,yiT,h,y^2T,K,^nT)■ 
n=l •' 

(87) 

Representing the exclusive particle density as the square of multiparticle amplitudes, 
expressing the configuration space amplitudes as Fourier transforms of amplitudes in 
momentum space, and performing integration over d'^X2T ■ ■ ■ d'^XnT finally one arrives 
at Eq. (^. 

13 Appendix H 

For our purposes let us quote the following identities for Bessel functions (|T^ and 
references therein) : 

oo 

MxOk'k) = E e''"^''-''"^Mxekp)Mx9,,p), (88) 



m=— oo 



oo 



e 



ikx cosip \ ^ -m imip 



E t'^e^'^^JUkx), (89) 



m=— oo 

oo 







, 6(a — a') 

dx xJq[x a)jQ[x a ) = ; (90) 



where 6k'k denotes the relative angle between momenta k' k, which in the reference 
frame with the z-axis placed along P direction takes form (^), and ipk'k = f — f', 
where the angles ip, (p' denote the azimuthal angles of vectors kx and k^ on the 
transversal plane respectively. 

Now let us repeat the scheme Eq.(^)— Eq.(0) iteratively in (^31). Then, after 
introducing the explicit form of Fourier transform we arrive at the expression : 

Pp (A;, xt(p),0)= (91) 
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(27r)2 

°° (26)" , „ (P\ 2b 1 



X y In" ^ n Qp{k)Qp{k') 

n\ \k) 27ck^epkepk' ^ ^ 

dXiXiJoiXiOk'k)-- dXnXnMXnOk'k)] 

x< / rfaiJo(xiai) da^Jo{xiar).. I danJo{xnan) L da^Jo{xnaJ > 

K. k k k k ) 

It looks still quite complicated, however one can notice that the term in brackets : 
1 / Qa daiJoi^xiai) I da^Jo{xiav).. L danJoixnttn) L ' daj^Jo{xna.^) > ; 



k k 



(92) 

which contains integration of 2n Bessel functions, simplifies significantly if one adds 
to (|9lD some extra off-diagonal terms, which are equal to in the quasi-diagonal 
limit. Namely : 

DLA 1 r^pk r^pk' I rSpk r^py ^ ^ 

{} = — daiJo{xiai) da^jQ^Xiay). . I danJoixnan)L da'^Jo{xnaJ. 

k k k k 

(93) 

Then expression (^) takes the simpler form : 

Pp (A:, xt(p),0) = ^ke{^ < A: < P)^ X 

X L dOkP L ddk'P / d^ / d^'e-^'^^kpXT CO, ^+ikeypXT cos ^ /g4^ 

JQ^ JQo Jq Jq 



k k 



5Z TT;2" ( T ) 1 / dxxJo{x6k'k) L daJoixa) I da'Jo{xa)> . 

Let us perform the integration over the angles (p, (p'. Replacing all the terms of the 
type JoixOk'k) by their expansions (0), expanding then Fourier transforms in terms 



of (|89D , and performing the integration over (p, ip' explicitly, we arrive at (|6ll) 



However, the part of expression (^) containing Bessel functions can be rewritten 



as 



n oo 



U $Z JmX^i^Pk)Jmi{xSpk') Jmi+. .+m„{xTkOpk) Jmi+. .+mn {XxkOpk' 
. j=l mi=—co 
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1 r2. 

— / difk'k 



27r Jo 



k'kj 



\i=l 



(95) 



^pfc' + ^"pk " 2 6'pfc/6'pfc cos(v9fc/fc), as usual. Substituting (p5|) into (|6T|) finally 
we obtain (158|). 



where yj^'A: denotes the relative angle between the vectors kj^, and 

(^k'k = 



14 Appendix I 

Let us investigate normalization of single particle density ( ^8]) in real space, i. e. the 
integral / (i^XTPp (fc, xt, 0). It reads : 

= 2bkQ{— < k < P) N ^ d6pk J^^ ^^"'^ dVk'kj dxTXTJoikOk^kXr) 



X Jo(^ 86 ln(— ) |y dxxJo{x9kk') Jq^ daJo(xa) da' Jo{xa') j) . (96) 

The term Jo{k9k'kXT) is the only part of expression (|96|) depending on xt- Integral 
dxTXTJo{k9k'kXT) can be performed, using identities (|88|) , (pO|) . It equals : 

dxTXTMkOk'kXT) = 27r5(v9fc/fc)^^%^— (97) 
k-^OkP 



Substituting ( pTl) into (|96D , and integrating over 6k' p and (pk'k, one arrives at : 

^b^.Qo , „^ dOpk 



J d'xrpT{k,^T,0) = J^(^<^<P)U 
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X 



Pk 



x/o(. 



A 86 1n(— )|y dxx daJo{xa) da'Jo{xa')>). (98) 

Performing the integration over x in the argument of function Iq with the help of 
identity (pO|), finally one obtains : 



dVp?(*,x.,0) ^ f e(^^ < * < P)/; |^/„(^8. ln(f (99) 
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which equals the normahzation of particle density in momentum space defined in 
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Figure captions 



Fig. 1 Feynman diagram for the production of m gluons in DLA, where P denotes 
the 4-momentum of the initial q{q) and ki denotes the 4-momentum of the ith pro- 
duced gluon. 



Fig. 2 Generating functional ([T9|) as a diagram series 



Fig. 3 Redefined kinematical regions in DLA. 

Fig. 4 Interference between different diagrams for d']¥{k[, ki). Remark: diagrams 
(a) and (b) are identical except of the position of ki {k[) leg. 

Fig. 5 Master equation for generating functional (p6D represented as a diagram 
series (for details see [[l^). Function -Pi',...,n';i,...,n introduces the parallel angular 
ordering (AO). 



Fig. 6 Function pp(k) from Eq. (|52D vs kx =\ ky | in Qq units for parameters 
b = 0.25, P/Qo = 243, = 1, k/Qo = 128, chosen following Pf. The power fit 



reads p'^ik) = (g;)"^'^^ * 0.00065 Qq - Plot (a) with the logarithmic scale for the 
vertical axis, plot (b) with the logarithmic scale for both vertical and horizontal axes. 

Fig. 7 Function p^p'^^\k,:Kt,0) from Eq. ( ^5]) vs. xt =\ | for parameters as in 



Fig. 6. The power fit reads pf{k, xt, 0) = {QoXt)'^-^^ * 4:.0E - 05 Qo- Plot (a) with 
the logarithmic scale for the vertical axis, plot (b) with the logarithmic scale for 
both vertical and horizontal axes. 

Fig. 8 Function pp(A;,Xf,0) from Eq. (p^) vs. Xt =| xt | for parameters as in 
Fig. 6. Plot (a) with the logarithmic scale for the vertical axis, plot (b) with the 
logarithmic scale for both vertical and horizontal axes. 
The power fit reads p'^(fc,xr,0) = (QoXt)'^-'^^ * 0.00018 Qo- 

Fig. 9 Diagramatic representation of the term 

/P\ ( f°° f^pk r^py 



26 In ^— j |y dxxJo{x6kk') J^^ daJo{xa) j^^ da'jQ^xa) 
dh — Assik,'k'T;k, kr); 
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The integration region lies in the overlap of conus Qs{k) and conus Qs{k')- 

Fig. 10 Diagramatic representation of p^p{k, xt, 0) from Eq. (^) as the chain of 
independent emission factors ( ]60| ) transformed onto the transverse xt plane by the 
Bessel factor of primary emission from the parent P, namely Jo(| — k'^ \ ^t)- 
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